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Abstract 

O . 

, Twenty five years have passed by since models of accretions and jets have separately 

■ emerged. Today, it is understood that these two objects are related to each other in 

a fundamental way. In a binary system, matter from an accretion disk enters into a 

On : black hole. A part of it is bounced back because of the centrifugal barrier, radiation 

pressure or magnetohydrodynamic effects, to form jets and bipolar outflows which carry 
away excess angular momentum. In the case of AGNs containing black holes, accretion 
disks form out of stellar winds and similar processes as above form cosmic radio jets. 
We present a general review of the study of the accretion disks and outflows in a 
I , coherent manner, especially emphasizing global inflow-outflow solutions (GIOS). We 

■ also present a few observational consequences of wind production from the accretion 
t-H ' disks on spectral properties of the accretion disks. 

00 '. 

\ To appear in the proceedings of the Mini Workshop on Fluid Dynamics at the Dacca Uni- 

versity (August, 1998). Ed. A. Hosain. 

Of 

O ■ 1 Introduction 

r/3 . Study of accretion disks and outflows around black holes began twenty- five years ago [1-2]. 

The subject has evolved considerably since then and it is now clear that these two apparently 
J> . dissimilar objects are related to each other. The accretion solutions of purely rotating disk 

' have been improved to include the effect of radiation pressure [3-7] (See Chakrabarti [8] 

• for a review.). Jet solutions have changed from speculative ideas such as de-Laval nozzles 

' [9] to electrodynamically acceleration model [10], self-similar centrifugally driven outflows 

[11], 'cauldrons' [12] etc. Centrifugally driven outflows are subsequently modified to include 
accretion disks [13]. Chakrabarti & Bhaskaran [14] (see also, Contopoulos, [15]) showed that 
it is easier produce outflows from a sub-Keplerian inflow. 

Parallelly, efforts were on to study the accretion and jets within the same framework. 
Chakrabarti [16-17] found that the same solution of purely rotating flows around a black 
hole could describe accretion flows on the equatorial plane and pre-jet matters near the 
axis. With a natural angular momentum distribution of l(r) = c\(r) n , (where c and n 
are constants and A is the von Zeipcl parameter) it was found that for large c and small 
n (n < 1), solutions are regular on the equatorial plane and they describe thick accretion 
disks. For small c and large n (n > 1), the solutions are regular on the axis and they 
describe pre-jet matters. It was speculated that some viscous process might be responsible 
to change the parameters from one set to the other. Even in the absence of viscosity, 
constant angular momentum flow was found to bounce back from the centrifugal barrier 
in numerical simulations of accretion flows [18]. Eggum, Coroniti & Katz [19] considered 
radiatively outflows emerging from a Keplerian disk. 

Further progress of this topic required a fundamental understanding of these flows which 
emerged in the late 1980s. It was found out just as Bondi solution [20] of accretion onto stars 
is related to outflow solutions [21], fundamentally transonic black hole accretions and winds 
are also related to each other. All possible accretion and wind type solutions are found, 
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including solutions which may contain standing shocks and the entire parameter space is 
classified according to the nature of solutions [22-23]. Fig. 1 (taken from Chakrabarti, [24]) 
shows these solutions (Mach number is plotted against logarithmic radial distance [in Units 
of the Schwarzschild Radius] and outer boxes, and specific energy is plotted against specific 
angular momentum in the central box) and the classification of the parameter space when 
the Kerr parameter a — 0.5 and when the equatorial plane solutions are considered (similar 
solutions are present for conical flows in winds as well, see [23]). The inward pointing arrows 
indicate accretion solutions and the outward pointing arrows indicate wind solutions. The 
flow from regions I, O pass through the inner sonic point and outer sonic point respectively. 
Those from NSA have no shock in accretion, from SA have shocks in accretion, from NSW 
have no shocks in winds, and SW have shocks in winds respectively. The Global Inflow 
Outflow Solutions (GIOS) as will be described in §3 combine one solution of each kind to 
produce wind out of accretion. The horizontal line in the central box denotes the rest mass of 
the inflow. Note that the outflows are produced only when the specific energy is higher than 
the rest mass energy. Flow with with lesser energy produces solutions with closed topologies 
(I* and O*). When viscosity was added the nature of the solutions changed fundamentally 
[25] (see, [26-27] for details) allowing matter to directly come out of a Kcplerian disk and 
enter into a black hole through the sonic point. These solutions, both for the accretion 
and winds have been verified by complete time-dependent simulations [28-31]. In the case 
of steady state solutions, outflows are found to occur between the centrifugal barrier and 
the funnel wall (Fig. 2a below) while in a non-steady solution the outflow could spread in 
regions outside the centrifugal barrier as well. 

2 Recent Progresses in Accretion Disk Solutions 

Fig. 2a shows the general picture that emerges of the non-magnetized inflow and outflow 
around a black hole. Centrifugal force tends to fight against gravity close to a black hole 
to produce a denser region called CENBOL (centrifugal barrier supported boundary layer). 
Matter farther out rotates in a Keplerian disk, but close to the black hole the flow is puffed 
up due to radiation pressure or ion pressure depending on whether the accretion rate is 
sufficiently high or not. The two-dimensional nature of the density distribution is given in 
Chakrabarti [24]. 

Chakrabarti & Titarchuk [32] using this generalized accretion solution suggested that 
hard and soft states of a black hole could be understood simply by re-distribution of matter 
between Keplerian and sub-Keplerian components which is effected by variation of viscosity 
in the flow in the vertical direction. This general conclusions is verified observationally 
(Zhang et al. 1997). The constancy of energy spectral index a e (F v oc v~ a ') separately 
in hard and soft states, as well as possible shifting of the inner edge with the Keplerian 
component with accretion rate as predicted by the general accretion solution [26] have also 
been verified [34-36] . It was also recognized that the non-stationary solutions of the accretion 
flows [37-38] might cause quasi-periodic oscillations. It was pointed out that the outflows, 
and not accretion disks, could be responsible for the Iron K a lines and the so-called reflection 
components [32]. A schematic representation of the detailed picture as above is drawn in 
Fig. 2b. We shall use this picture below to analytically estimate the outflow rates. 

A new concept which was found useful in identifying black holes by using spectral features 
alone is the 'bulk motion Comptonization'. Basically, as matter flows into a black hole 
rapidly with almost the speed of light, photons scattering from them gain momentum and 
frequency is blue-shifted due to Doppler effect [32, 39]. The resulting spectra is of power- 
law and extends till about an MeV. The presence of this part of the spectra in soft states 
is widely accepted to be the only way to identify a black hole most convincingly, whether 
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galactic or extragalactic, since the all absorbing property of the horizon is directly used in 
obtaining the spectral slope. 

Black holes in X-ray binaries often show quiescence states. This could be the result of 
very low accretion rates. A novel solution that was proposed by Das & Chakrabarti [40] (also 
see, [41-42]) is that outflows generated from the inflow could, in some certain circumstances 
especially when the accretion rate is low, be so high that the disk may be almost evacuated. 
They proposed that such outflows could generate quiescence states of a black hole candidate. 

3 Global Inflow- Outflow Solutions (GIOS) 

Outflows are common in many astrophysical systems which contain black holes and neu- 
tron stars. Difference between stellar outflows and outflows from these systems is that the 
outflows in these systems have to form out of the inflowing material only. We now present 
a simple analytical approach by which the outflow rate is computed out of the inflow rate. 
The ratio of these two rates is found to be a function the compression ratio of the gas at 
the boundary between CENBOL and the inflow. 

The problem of mass outflow rate in the context of a black hole accretion has been 
attacked quite recently [27, 43]. A simple approach widely used in stellar astrophysics has 
been followed where the flow upto the sonic point is assumed to be isothermal. This was 
possible due to the novel understanding that around a black hole, the centrifugal pressure 
supported dense matter could behave like a 'boundary layer'. This CENBOL is hot, puffed 
up and very similar to the classical thick accretion disk, except that at the inner edge 
the matter is strongly advective. The thermal pressure is expected to drive matter out in 
between the centrifugal barrier and the funnel wall [31] as in Fig 2a. However, we use a 
simple model (Fig. 2b) where both the inflow is axially symmetric and wedge shaped while 
the outflow is conical with solid angles 6j„ and Q ou t respectively. 

CENBOL could form either with or without shocks as long as the angular momentum 
of the flow close to the compact object is roughly constant as is the case in reality [22-24]. 
This region replaces [32, 44-45] the so called 'Compton cloud' in explaining hard and soft 
states of black hole. The oscillation of this region can explain the general properties of the 
quasi-periodic oscillation [37,38,46] from black holes and neutron stars. It is therefore of 
curiosity if this region plays any major role in formation of outflows. 

Several authors have also mentioned of denser regions closer to a black hole due to differ- 
ent physical effects. Chang & Ostriker [47] showed that pre-heating of the gas could produce 
standing shocks at a large distance (typically a few tens of thousands Schwarzschild radii 
away). Kazanas & Ellison [48] mentioned that pressure due to pair plasma could produce 
standing shocks at smaller distances around a black hole as well. Our computation is insensi- 
tive to the actual mechanism by which the boundary layer is produced. All we require is that 
the gas should be hot at the region where the compression takes place. Thus, since Comp- 
tonization processes cool this region [32] for larger accretion rates (M £ 0-^MEddington) our 
process is valid only for low-luminosity objects, consistent with current observations. Begel- 
man & Rees [12] talked about a so-called 'cauldron' model of compact objects where jets 
were assumed to emerge from a dense mixture of matter and radiation by boring a de-Laval 
nozzle as in Blandford & Rees [9] model. The difference between this model and the present 
one is that very high accretion rate was required {Mm ~ 1000Mb) there while we con- 
sider thermally driven outflows out of accretion with smaller rates. Second, the size of the 
'cauldron' was thousands of Schwarzschild radii (where gravity was so weak that channel 
has to have shape of a de-Laval nozzle to produce transonicity), while we have a CENBOL 
of about 10Rg (where the gravity plays an active role in creating the transonic wind) in 
our mind. Third, in the present case, matter is assumed to pass through a sonic point 
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using the pre-determined funnel where rotating pre-jet matter is accelerated [16] and not 
through a 'bored nozzle' even though symbolically a quasi-spherical CENBOL is considered 
for mathematical convenience. 

Once the presence of the centrifugal pressure supported boundary layer (CENBOL) is 
accepted, the mechanism of the formation of the outflow becomes clearer. One basic criteria 
is that the outflowing winds should have positive Bernoulli constant [22]. Just as photons 
from the stellar surface deposit momentum on the outflowing wind and keeps the flow 
roughly isothermal [49] at least upto the sonic point, one may assume that the outflowing 
wind close to the black hole is kept isothermal due to deposition of momentum from hard 
photons. In the case of the sun, it's luminosity is only 10~ 5 LEdd and the typical mass 
outflow rate from the solar surface is 1O _14 M year -1 [50]. Proportionately, for a star with 
a Eddington luminosity, the outflow rate would be 10 _9 M© year -1 . This is roughly half as 
the Eddington rate for a stellar mass star. Thus if the flow is compressed and heated at 
the centrifugal barrier around a black hole, it would also radiate enough to keep the flow 
isothermal (at least up to the sonic point) if the efficiency were exactly identical. Physically, 
both requirements may be equally difficult to meet, but in reality with photons shining 
on outflows near a black hole with almost An solid angle (from funnel wall) it is easier to 
maintain the isothermality in the slowly moving (subsonic) region in the present context. 
Another reason is this: the process of momentum deposition on electrons is more efficient 
near a black hole. The electron density n e falls off as r -3 / 2 while the photon density n 7 falls 
off as r~ 2 . Thus the ratio n e /n 7 cx r 1 / 2 increases with the size of the region. Thus a compact 
object will have lesser number of electrons per photon and the momentum transfer is more 
efficient. In a simpler minded way the physics is scale- invariant, though. In solar physics, it 
is customary to chose a momentum deposition term which keeps the flow isothermal to be 
of the form [51], 

F r = I Ddr 

J R s 

where, D is the momentum deposition (localized around r p ) factor with a typical spatial 
dependence, 

D = D e- a(r/r ?- 1)2 

Here, Dq, a are constants and R s is the location of the stellar surface. Since r and r p comes 
in ratio, exactly same physical consideration would be applicable to black hole physics, 
with the same result provided Dq is scaled with luminosity (However, as we showed above, 
Do goes up for a compact object due to higher solid angle.). However, as Chakrabarti & 
Titarchuk [32] showed, high accretion rate (M ^ 0.3MEdd ) will reduce the temperature of 
the CENBOL catastrophically, and therefore our assumption of isothermality of the outflow 
would severely breakdown at these high rates simply because cooler outflow would have 
sonic point very far away from the black hole. It is to be noted that in the context of stellar 
physics, it is shown [52] that the temperature stratification is the outflowing wind has little 
effect on the mass loss rate. Effect of radiation momentum deposition on the topology of 
the outflows is separately discussed in Chattapadhyay [53] . 

3.1 Derivation of the outflow rate using simple GIOS 

The accretion rate of the incoming accretion flow is given by, 

Min = ^npVr 2 . (1) 

Here, 0, n is the solid angle subtended by the inflow, p and d are the density and velocity 
respectively, and r is the radial distance. For simplicity, we assume geometric units (G = 
1 = Mbh = c; G is the gravitational constant, Mbh is the mass of the central black hole, 
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and c is the velocity of light) to measure all the quantities. In this unit, for a freely falling 
gas, 

Hr) = l 1 -^} 112 (2) 

and 



p(r) = ^(l-T)-^r-^ (3) 

Here, T/r 2 (with T assumed to be a constant) is the outward force due to radiation. 

We assume that the boundary of the denser cloud is at r = r s (typically a few Schwarzschild 
radii, see, Chakrabarti [24]) where the inflow gas is compressed. The compression could be 
abrupt due to standing shock or gradual as in a shock-free flow with angular momentum. 
This details are irrelevant. At this barrier, then 

P+{r s ) = Rp-(r s ) (4a) 

and 

^ + (r s ) = i?" 1 1 ?_(r s ) (46) 

where, R is the compression ratio. Exact value of the compression ratio is a function of the 
flow parameters, such as the specific energy and the angular momentum (e.g., [22-24]) Here, 
the subscripts — and + denote the pre-shock and post-shock quantities respectively. At the 
shock surface, the total pressure (thermal pressure plus ram pressure) is balanced. 

P-{r s ) + p-{r s )$ 2 _{r s )=P + {r s )+ P+ {r s )$l{r s ). (5) 

Assuming that the thermal pressure of the pre-shock incoming flow is negligible compared 
to the ram pressure, using eqs. 4(a-b) we find, 

P + (r s ) = ^p-(r s )0 2 _(r s ). (6) 
The isothermal sound speed in the post-shock region is then, 

2 p+ (fl-i)(i-r) i (i-r) 

s P+ R 2 r s f r s [i> 

where, /o = R 2 /(R — 1). An outflow which is generated from this dense region with very 
low flow velocity along the axis is necessarily subsonic in this region, however, at a large 
distance, the outflow velocity is expected to be much higher compared to the sound speed, 
and therefore the flow must be supersonic. In the subsonic region of the outflow, the 
pressure and density are expected to be almost constant and thus it is customary to assume 
isothermality condition up to the sonic point [49]. With isothcrmality assumption or a 
given temperature distribution (T oc r~P with [3 a constant; see eq. [22] below) the result 
is derivable in analytical form. The sonic point conditions are obtained from the radial 
momentum equation, 

d§ ldP i-r „ 
^— + -— + ^=0 8 

dr p dr r A 



and the continuity equation 



r 2 dr 

in the usual way, i.e., by eliminating dp/dr, 

dti N 



1 d(p#r 2 ) = Q 



dr D ^ 
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where 

Ar 2Cl 1 - r 
N=—?-- — — 

and 

r 2 

and putting N = and D = conditions. These conditions yield, at the sonic point r = r c , 
for an isothermal flow, 

0(r c ) = C a . (11a) 

and 

where, we have utilized eq. (7) to substitute for C s . 

Since the sonic point of a hotter outflow is located closer to the black hole, clearly, the 
condition of isothermality is best maintained if the temperature is high enough. However if 
the temperature is too high, so that r c < r s , one has to solve this case more carefully, using 
considerations of Fig. 2a, rather than of Fig. 2b. 

The constancy of the integral of the radial momentum equation (eq. 8) in an isothermal 
flow gives: 

C 2 s ln p + - ifl = \d + C 2 Jn Pc -— (12) 

where, we have ignored the initial value of the outflowing radial velocity i9(r s ) at the dense 
region boundary (r — r s ), and also used eq. (11a). We have also put p(r c ) = p c and 
p(r a ) = p+. Upon simplification, we obtain, 

p c = p+exp(-f) (13) 

where, 

/ = fo - 2 

Thus, the outflow rate is given by, 

M out = Q out p c $ c rl (14) 

where, <d ou t is the solid angle subtended by the outflowing cone. Upon substitution, one 
obtains, 

%*=B* = JT i jf S o /a exp{-f) (15) 
which, explicitly depends only on the compression ratio: 

~k~ = ^ = e-I [ ^T ] exp{ 2 - r—i ] (16) 

apart from the geometric factors. Notice that this simple result does not depend on the 
location of the sonic points or the the size of the dense cloud or the outward radiation force 
constant T. This is because the Newtonian potential was used throughout and the radiation 
force was also assumed to be very simple minded (T/r 2 ). Also, effects of centrifugal force 
was ignored. Similarly, the ratio is independent of the mass accretion rate which should 
be valid only for low luminosity objects. For high luminosity flows, Comptonization would 
cool the dense region completely [32] and the mass loss will be negligible. Pair plasma 
supported quasi-spherical shocks forms for low luminosity as well [48]. In reality there 
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would be a dependence on these quantities when full general relativistic considerations of 
the rotating flows are made. Exact and detailed computations using both the transonic 
inflow and outflow (where the compression ratio R is also computed sclf-consistently) are 
presented elsewhere [40] . 

Figure 3 contains the basic results. The solid curve shows the ratio R m as a function 
of the compression ratio R (plotted from 1 to 7), while the dashed curve shows the same 
quantity as a function of the polytropic constant n = (7 — I) -1 (drawn from n = 3/2 to 3), 7 
being the adiabatic index. The solid curve is drawn for any generic compression ratio and the 
dashed curve is drawn assuming the strong shock limit only: R — (7 + 1 ) / (7 — I ) = 2n + 1 . In 
both the curves, O ut ~ @m has been assumed for simplicity. Note that if the compression 
does not take place (namely, if the denser region does not form), then there is no outflow 
in this model. Indeed for, R = 1, the ratio Rm is zero as expected. Thus the driving force 
of the outflow is primarily coming from the hot and compressed region. The basic form 
is found to agree with results obtained from rigorous calculation with transonic inflow and 
outflow. 

In a relativistic inflow or for a radiation dominated inflow, n — 3 and 7 = 4/3. In the 
strong shock limit, the compression ratio is R = 7 and the ratio of inflow and outflow rates 
becomes, 

R m = 0.052 (17a) 

For the inflow of a mono-atomic ionized gas n = 3/2 and 7 = 5/3. The compression ratio is 
R = 4, and the ratio in this case becomes, 

Rm = 0.266 (176) 

Since /o is smaller for 7 = 5/3 case, the density at the sonic point in the outflow is much 
higher (due to exponential dependence of density on /o, see, eq. 7) which causes the higher 
outflow rate, even when the actual jump in density in the postshock region, the location of 
the sonic point and the velocity of the flow at the sonic point are much lower. It is to be noted 
that generally for 7 > 1.5 shocks are not expected [23], but the centrifugal barrier supported 
dense region would still exist. As is clear, the entire behavior of the outflow depends only 
on the compression ratio, R and the collimating property of the outflow otl t/@m- 

Outflows are usually concentrated near the axis, while the inflow is near the equatorial 
plane. Assuming a half angle of 10° in each case, we obtain, 

e - 2 < e 

w in g ) ^out 



and 



B ""' * (18) 



& in 36 

The ratios of the rates for 7 = 4/3 and 7 = 5/3 are then 

R m = 0.0045 (19a) 

and 

R m = 0.023 (196) 

respectively. Thus, in quasi-spherical systems, in the case of strong shock limit, the outflow 
rate is at the most a couple of percent of the inflow. If this assumption is dropped, then for 
a cold inflow, the rate could be higher by about fifty percent (see, Fig. 3). 

It is to be noted that the above expression for the outflow rate is strictly valid if the 
flow could be kept isothermal at least up to the sonic point. In the event this assumption 
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is dropped the expression for the outflow rate becomes dependent on several parameters. 
As an example, we consider a polytropic outflow of same index 7 but of a different entropy 
function K (We assume the equation of state to be P — Kp 7 , with 7 7^ 1). The expression 
(lib) would be replaced by, 

for s 



C 2 7 

and eq. (12) would be replaced by, 



(20) 



i-r ,1 N , l-r 

na 2 + - — = (-+n)a 2 - — (21) 

where n = 1/(7 — 1) is the polytropic constant of the flow and a + = (■jP + /p + ) 1 ^ 2 and 
a c = {'j P c / Pc) 1 ^ 2 are the adiabatic sound speeds at the starting point and the sonic point of 
the outflow. It is easily shown that a power law temperature fall off of the outflow (T cx r"' 3 ) 
would yield 

RrH = ^&nh^, (22) 

where, Ki and K Q are the entropy functions of the inflow and the outflow. This derivation 
is strictly valid for a non-isothermal flow. Since Ki < K Q , n > 3/2 and /o, for Q ou t ~ &i n , 
R,n « 1 is guaranteed provided (3 > |, i.e., if the temperature falls for sufficiently rapidly. 
For an isothermal flow (3 = and the rate tends to be higher. Note that since n ~ 00 in this 
case, any small jump in entropy due to compression will off-balance the the effect of f$ 3 ^ 2 
factor. Thus Rm remains smaller than unity. The first factor decreases with entropy jump 
while the second factor increases with the compression ratio (R) when (3 < 3/2. Thus the 
solution is still expected to be similar to what is shown in Fig. 3. 



3.2 Results from Rigorously obtained GIOS 

When inflow and outflow solutions are obtained more rigorously, i.e., actually making the 
flow pass through separate sonic surfaces, one has to include the effect of mass-loss on the 
Rankinc-Hugoniot condition at the boundary between CENBOL and the accretion flow [40]. 
Accordingly, we use, 

M+ = (1 - Rin)M- (23) 

where, the subscripts + and — denote the pre- and post-shock values respectively. Since due 
to the loss of matter in the post-shock region, the post-shock pressure goes down, the shock 
recedes backward for the same value of incoming energy, angular momentum & polytropic 
index. The combination of three changes, namely, the increase in the cross-sectional area 
of the outflow and the launching velocity of the outflow and the decrease in the post-shock 
density decides whether the net outflow rate would increased or decreased than from the 
case when the exact Rankinc-Hugoniot relation was used. 

Fig. 4 shows a typical global inflow-outflow solution (GIOS) where actual transonic prob- 
lem was solved. The input parameters arc £ = 0.0005, A = 1.75 and 7 = 4/3 corresponding 
to relativistic inflow. The solid curve with an arrow represents the pre-shock region of the 
inflow and the long-dashed curve represents the post-shock inflow which enters the black 
hole after passing through the inner sonic point (I). The solid vertical line at X s ^ (the left- 
most vertical transition; the notation of the shock location is from [22]) with double arrow 
represents the shock transition obtained with exact Rankine-Hugoniot condition (i.e., with 
no mass loss). The actual shock location obtained with modified Rankinc-Hugoniot condi- 
tion (eq. 23) is farther out from the original location X s3 . Three vertical lines connected 
with the corresponding dotted curves represent three outflow solutions for the parameters 
7o = 1.3 (top), 1.15 (middle) and 1.05 (bottom). The outflow branches shown pass through 
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the corresponding sonic points. It is evident from the figure that the outflow moves along 
solution curves which are completely different from that of the 'wind solution' of the inflow 
which passes through the outer sonic point '0\ The mass loss ratio R m in these cases are 
0.256, 0.159 and 0.085 respectively. This figure is taken from [40]. Fig. 3 of that work 
shows the general behaviour of the variation of with compression ratio of the shock. 
This agrees well with what is obtained analytically (Fig. 3 of this paper) . 

Das & Chakrabarti [40], for the first time, pointed out that in certain region of the 
parameter space, the mass loss could be so much that disk would be almost evacuated and 
conjectured that inactive phases of a black hole (such as the present day Sgr A* at our 
galactic center) may be formed when such profuse mass loss takes place. 

It is to be noted that although the existence of outflows are well known, their rates are 
not. The only definite candidate whose outflow rate is known with any certainty is probably 
SS433 whose mass outflow rate was estimated to be M ou t ~ 1-6 x 10~ 6 f~ 1 n^ 3 D^Mq yr _1 
[54], where / is the volume filling factor, nis is the electron density n e in units of 10 13 
cm" 3 , Dc, is the distance of SS433 in units of 5kpc. Considering a central black hole of 
mass 10M Q , the Eddington rate is Me<i ~ 0.2 x 10~ 7 M Q yr _1 and assuming an efficiency 
of conversion of rest mass into gravitational energy r\ ~ 0.06, the critical rate would be 
roughly M cr i t = MEd/v ~ 3.2 x 10~ 7 M Q yr _1 . Thus, in order to produce the outflow rate 
mentioned above even with our highest possible estimated Rm ~ 0.4 (see, Fig. 2), one must 
have Mi n ~ 12.5M crit which is very high indeed. One possible reason why the above rate 
might have been over-estimated would be that below 10 12 cm from the central mass [54], 
^13 >> 1 because of the existence of the dense region at the base of the outflow. 

In numerical simulations the ratio of the outflow and inflow has been computed in several 
occasions [19,29]. Eggum et al. [19] found the ratio to be R m ~ 0.004 for a radiation 
pressure dominated flow. This is generally comparable with what we found above (eq. 19a). 
In Molteni et al. [29] the centrifugally driven outflowing wind generated a ratio of Rm ~ 0.1. 
Here, the angular momentum was present in both inflow as well as outflow, and the shock 
was not very strong. Thus, the result is again comparable with what we find here. 

3.3 Spectral Softening in Presence of Profuse Mass Loss 

Creation of winds out of accretion disk has an interesting consequence. As more fraction 
of matter is expelled, it becomes easier for the soft photons from the Keplerian disk (see 
Fig. 2a) surrounding the advective region to cool this region due to Comptonization (see, 
Chakrabarti & Titarchuk for details). In other words, the presence of winds would soften 
the spectra of the power-law component for the same multicolour blackbody component. 
Such an observation would point to profuse mass loss from the hot advective region. 

The model we use here is the two component accretion flow which has a CENBOL 
close to the hole. We assume a weakly viscous flow of constant specific angular momentum 
A = 1.75 (in units of 2GM/c) for the sake of concreteness. The Keplerian component close 
to the equatorial plane has a low accretion rate (M,„ ~ 0.05 — 0.3 in units of the Eddington 
rate) and the sub-Keplerian halo surrounding it has a higher rate (Mh ~ 1 in units of 
Eddington rate). Before the accreting matter hits the inner advective region, both the rates 
are constant, but as Das & Chakrabarti [40] has shown, winds, produced from CENBOL 
will deplete the disk matter at the rate determined by the temperature of the CENBOL, 
when other parameters, such as the specific angular momentum and specific energy are kept 
fixed. 

Figures 5 and 6 show the outcome of our calculation of the spectra for three different 
accretion rate of the Keplerian component Mj„ [43] . The mass of the central black hole is 
chosen to be M = 10M©. The size of the CENBOL is assumed to be 10r s (where r g is 
the Schwarzschild radius), a typical location for the sub-Keplerian flow of average specific 
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angular momentum A = 1.75 and specific energy £ = 0.003. Following Das & Chakrabarti 
[40], we first compute the mass outflow rate from the advective region. The long dashed 
curve in Fig. 5 shows the variation of the percentage of mass loss (vertical axis on the right) 
as a function of the inflow accretion rate. The dotted curve and the solid curve denote the 
variation of the energy spectral index a (F u oc u~ a ) with and without winds taken into 
account. Note that the spectra is overall softened (a increased) when winds are present. 
For higher Keplerian rates, the mass loss through winds is negligible and therefore there is 
virtually no change in the spectral index. For lower inflow rates, on the other hand, mass 
loss rate is more than twenty percent. It is easier to Comptonize the depleted matter by 
the same number of incoming soft photons and therefore the spectra is softened. 

In Fig. 6, we show the resulting spectral change. As in Fig. 5, solid curves represent 
solutions without winds and the dotted curves represent solutions with winds. Solid curves 
are drawn for M,„ = 0.3 (uppermost at the bump), 0.15 (middle at the bump) and 0.07 
(lowermost at the bump) respectively. For Mi n = 0.3 both curves are identical. Note the 
crossing of the solid curves at around 10 1S 6 Hz (15 keV) when winds are absent. This is 
regularly observed in black hole candidates. If this is shifted to higher energies, the presence 
of winds may be indicated. 

Strong winds are suspected to be present in Sgr A* at our Galactic Center (see, Genzel 
et al. [55] for a review, and Eckart & Genzel [56]). Chakrabarti [9] suggested that the 
inflow could be of almost constant energy transonic flow, so that the emission is inefficient. 
However, from global inflow-outflow solutions [GIOS], Das & Chakrabarti [40] showed that 
when the inflow rate itself is low (as is the case for Sgr A*; ~ 10~ 3 to 10~ 4 MEddington) the 
mass outflow rate is very high, almost to the point of evacuating the disk. This prompted 
them to speculate that spectral properties of our Galactic Center could be explained by 
inclusion of winds. This will be done in near future. Not only our Galactic Center, the 
consideration should be valid for all the black hole candidates (e.g., V404 Cyg) which are 
seen in quiescence. 

Chakrabarti & Titarchuk [32] suggested that the iron K a line as well as the so called 
'reflection component' could be due to outflows off the advective region. Combined with the 
present work, we may conclude that simultaneous enhancement of the 'reflection component' 
and/or iron K a line intensity with the softening of the spectra in hard X-rays would be a 
sure signature of the presence of significant winds in the advective region of the disk. 
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FIGURE CAPTIONS 



Fig. 1 : Classification of the parameter space (central box) in the energy- angular momen- 
tum plane in terms of various topology of the black hole accretion. Eight surrounding 
boxes show the solutions from each of the independent regions of the parameter space. 
Each small box shows Mach number M against the logarithmic radial distance r 
(measured in units of GM BH / c 2 ) Contours are of constant entropy accretion rate M. 
Similar classification is possible for all adiabatic index 7 < 1.5. For 7 > 1.5, only 
the inner sonic point is possible other than an unphysical 'O' type point (Chakrabarti 
1996c). 

Fig. 2a : Cartoon diagram of a very general inflow-outflow configuration of non-magnetized 
matter around a compact object. Keplerian and sub-Keplerian matter accretes along 
the equatorial plane. Centrifugally and thermally driven outflows preferentially emerge 
between the centrifugal barrier and the funnel wall. 

Fig. 2b Schematic diagram of inflow and outflow around a compact object. Hot, dense 
region around the object either due to centrifugal barrier or due to plasma pressure 
effect or pre-heating, acts like a 'stellar surface' from which the outflowing wind is 
developed. 

Fig. 3 Ratio Rm of the outflow rate and the inflow rate as a function of the compression 
ratio of the gas at the dense region boundary (solid curve). Also shown in dashed curve 
is its variation with the polytropic constant n in the strong shock limit. Solid angles 
subtended by the inflow and the outflow are assumed to be comparable (Chakrabarti 
1998a). 

Fig. 4 : Variation of the percentage of mass loss (long dashed curve and right axis) and 
the energy spectral index a (F v oc a~ a ) (solid and dotted curves and left axis) with 
the accretion rate (in units of Eddington rate) of the Keplerian component. Solid 
curve is drawn when winds are neglected from the advective region and dotted curve 
includes effect of winds. Overall spectra is softened when the inflow rate is reduced 
(Chakrabarti 1998b). 

Fig. 5 : Spectra of emitted radiation from the accretion disk with (dotted) and without 
(solid) effects of winds. Hard X-ray component is softened while keeping soft X-ray 
bump unchanged (Chakrabarti 1998b). 

Fig. 6 : Few typical solutions which combine accretion and outflow. Input parameters are 
£ = 0.0005, A = 1.75 and 7 = 4/3. Solid curve with an incoming arrow represents the 
pre-shock region of the inflow and the dashed curve with an incoming arrow represents 
post-shock inflow which enters the black hole after passing through the inner sonic 
point (I). Dotted curves are the outflows for various j (marked). Open circles are 
sonic points of the outflowing winds and the crossing point 'O' is the outer sonic 
point of the inflow. The leftmost shock transition (X s3 ) is obtained from unmodified 
Rankine-Hugoniot condition, while the other transitions are obtained when the mass- 
outflow is taken into account. 
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